Effect of interface states on spin-dependent tunneling in Fe/MgO/Fe tunnel junctions 
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The electronic structure and spin-dependent tunneling in epitaxial Fe/MgO/Fe(001) tunnel junc- 
tions are studied using first-principles calculations. For small MgO barrier thickness the minority- 
spin resonant bands at the two interfaces make a significant contribution to the tunneling con- 
ductance for the antiparallel magnetization, whereas these bands are, in practice, mismatched by 
disorder and/or small applied bias for the parallel magnetization. This explains the experimentally 
observed decrease in tunneling magnetoresistance (TMR) for thin MgO barriers. We predict that 
a monolayer of Ag epitaxially deposited at the interface between Fe and MgO suppresses tunneling 
through the interface band and may thus be used to enhance the TMR for thin barriers. 
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Magnetic tunnel junctions (MTJs) are miniature de- 
vices which consist of two ferromagnetic electrodes sepa- 
rated by an insulating barrier. These junctions are made 
in such a way that their magnetization may be switched 
between parallel and antiparallel states under the influ- 
ence of external magnetic field. This switching is ac- 
companied by an abrupt change of the electric conduc- 
tance of the MT J . MTJs aroused much attention due 
to their potential application in magnetic random-access 
memories and magnetic field sensors. In practical terms, 
the figure of merit is the tunneling magnetoresistance 
(TMR), which is defined by TMR = (G P - G A p)/G A p, 
where Gp and Gap are the conductances measured when 
the electrodes are magnetized parallel or antiparallel to 
each other. Recent reviews of spin-dependent tunneling 
in MTJs may be found in Refs. Ell 

Since the first observation of reproducible TMR 0, 
the majority of measurements were performed for amor- 
phous or polycrystalline barriers, most commonly AI2O3. 
The highest TMR values achieved for AI2O3 barriers were 
about 70% at room temperature 0. Meanwhile, the- 
oretical calculations based on layer KKR |f| and tight- 
binding methods Q predicted that much larger TMR val- 
ues may be obtained for coherent tunneling in epitaxial 
Fe/MgO/Fe(001) junctions due to strong spin filtering. 
The latter is enforced by the wave-function symmetry 
and its relation to the complex band structure of the 
barrier || • Very large TMR values exceeding 200% were 
indeed measured for such junctions by Parkin et al. 
and Yuasa et al. ^(J- Recently, a more accurate calcu- 
lation ^3 based on the FLAPW method confirmed the 
conclusions of Refs. 0,0. 

For device applications of MTJs it is critical to make 
the tunneling barrier as thin as possible in order to match 
the resistance of MTJs to other electronic components. 
Measurements for epitaxial Fe/MgO/Fe junctions show, 
however, that TMR decreases precipitously for barrier 
thickness below 2 nm [lfj. A detailed characterization 
of the MgO structure grown on Fe(001) single crystals 
demonstrates a pseudomorphic growth of MgO up to 6 
monolayers (ML) (« 1.2 nm), with misfit dislocations 
being formed for thicker films |l2j • The two latter exper- 



imental observations suggest that in the range of MgO 
thickness at which one might expect a ballistic tunnel- 
ing mechanism for conduction with no contribution from 
defect scattering, TMR drops down with decreasing the 
barrier thickness. The origin of this behavior is unknown. 
Also, these experimental facts are in disagreement with 
large values of TMR calculated for thin MgO barriers 



In this paper we demonstrate that the reduction of 
TMR in epitaxial Fe/MgO/Fe(001) junctions at small 
barrier thickness is controlled by the minority-spin inter- 
face band. The presence of this band was experimen- 
tally proven by Tiusan et al. 0]. We show that the 
transmission through this resonant channel is enhanced 
dramatically at small barrier thickness making a large 
contribution to the conductance in the antiparallel con- 
figuration and to the minority-spin conductance in the 
parallel configuration. The latter is, however, so sensi- 
tive to the mismatch in the potential at the two interfaces 
that it is, in practice, destroyed by disorder and/or ap- 
plied bias. This explains the sizable decrease in TMR 
for thin MgO barriers which is observed experimentally 
(lcj . We predict that a monolayer of Ag epitaxially de- 
posited at the interface between Fe and MgO suppresses 
tunneling through this interface band and may thus be 
used to enhance TMR for thin barriers. This provides a 
new way to make MTJs with a low resistance and high 
TMR that are required for device applications. In addi- 
tion, Ag interlayers protect the ferromagnetic electrodes 
from oxidation which is detrimental to TMR 0] . 

We calculate the electronic structure and tunneling 
conductance of Fe/MgO/Fe(001) MTJs with or without 
Ag interlayers using a tight-binding linear muffin-tin or- 
bital method (TB-LMTO) in the atomic sphere approx- 
imation (ASA) and the local density approximation 
(LDA) for the exchange-correlation energy. We use a full- 
potential LMTO (FP-LMTO) method to check the 
correctness of the ASA in describing the band structure 
of the MTJ. The principal-layer Green's function tech- 
nique is applied to calculate the conductance [17j . The 
atomic structure of the Fe/MgO/Fe junctions is taken 
from Ref. 0. To represent the electronic structure within 
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the ASA we use atomic spheres as described in Ref. 
The quality of this choice of the spheres is tested against 
our FP-LMTO calculations . In general, we find very 
good agreement between the ASA and FP results. In 
particular, the band offset between Fe and MgO is repro- 
duced very well. The Fermi level lies approximately 3.4 
eV above the MgO valence band maximum. 

The presence of the interface band can be visualized 
using the local density of states (DOS) resolved by trans- 
verse wave vector k||. Fig. ^ shows the minority-spin k||- 
resolved DOS at and around the Fermi energy, Ep, for 
the interface Fe layer in a Fe/MgO /Fe(OOl) junction. The 
band of states which is clearly seen in red in this figure is 
the interface band, absent in bulk Fe. This interface band 
can also be seen in the energy-resolved DOS as the nar- 
row peak near the Fermi level for minority-spin electrons 
(see, e.g., Fig. la in Ref. and Fig. 3b in Ref. @). 

It is known that properties of interface (surface) states 
depend on whether they are coupled to the bulk states or 
not [13] ■ In Fig. the interface states located about one 
quarter of the Brillouin zone width from the T point are 
interface resonances: They lie within the continuum of 
bulk Bloch states and therefore have a finite linewidth. 
On the other hand, the states forming two parallel curves 
in the corners of the Brillouin zone (inside the dark blue 
regions) are pure interface states. In the dark blue re- 
gions the DOS is zero in the bulk, and the interface 
states have zero linewidth. To resolve these states, we 
added an imaginary part of 10 -5 Ry to the energy. The 
two parallel bands correspond to bonding and antibond- 
ing combinations of the interface states localized at the 
two sides of the barrier [2l| . Near the points where these 
bands enter the bulk continuum and become resonances 
one can see strong peaks in the interface DOS, similar to 
those predicted within a simple tight-binding model • 

In a single-particle approximation the interface states 
projected into bulk band gaps do not contribute to the 
tunneling conductance. A possible way to include the 
contribution of such states was suggested by Ishida et 
al. 0. In our case, however, the interface resonances 
lie much closer to the T point compared to the pure in- 
terface states. Therefore, the resonances dominate the 
conductance, and the use of the single-particle approxi- 
mation does not lead to appreciable errors. 

A notable feature of the interface band at the Fe/MgO 
interface is its weak dispersion. This causes a significant 
change in the location of this band within the first Bril- 
louin zone when energy is shifted by a tiny amount of 
0.02 eV, as is seen in Fig.^ This feature makes any cal- 
culation of the interface states in Fe/MgO/Fe unreliable 
in terms of their Fermi level intercepts: The LDA itself 
does not provide 0.01 eV accuracy. It is very likely that 
this particular feature of the interface states is the rea- 
son why earlier calculations based on different methods 
6, 7, 11] result in very dissimilar shapes of the minority- 
spin conductance plotted as a function of k|| . 
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FIG. 1: Normalized minority-spin k||-resolved DOS at the 
interfacial Fe layer in Fe/MgO/Fe(001) MTJ for three values 
of energy: (a) 0.02 eV below E F , (b) at E F , and (c) 0.02 eV 
above E F . The scale is logarithmic. 



Panels (a)-(c) in Fig. [21 show the spin-resolved trans- 
mission for the Fe/MgO/Fe junction with 4 MLs of MgO 
for parallel and antiparallel magnetization. As is clearly 
seen from panel (b), the resonant interface band enhances 
the transmission in the minority spin channel. This en- 
hancement is most pronounced for small barrier thick- 
ness, because the interface band lies away from the T 
point, and therefore the resonant contribution to the 
transmission decays faster with barrier thickness com- 
pared to non-resonant. We find that for MgO thickness 
smaller than 6 MLs the contribution from minority-spin 
electrons in the parallel configuration becomes higher 
than that from majority-spin electrons. We note that 
in the calculation by Butler et al. Q this crossover does 
not occur down to 4 MLs of MgO, although the similar 
tendency is clearly seen from Fig. 16 in that paper. This 
disagreement likely results from the interface band cross- 
ing the Fermi level at a larger distance from the T point 
compared to our calculation. 

An important property of the minority-spin interface 
resonances is that they strongly contribute to the con- 
ductance in the parallel configuration only for ideal, sym- 
metric junctions, and only at zero bias. Indeed, it is seen 
in Fig.^that the interface DOS for these resonances ex- 
ceeds the DOS for neighboring regions of the surface Bril- 
louin zone by one to two orders of magnitude. Therefore, 
the interface resonances generate large tunneling current 
only if they match similar resonances at the other side 
of the barrier. As follows from Fig. ^ a bias voltage of 
the order of 0.01 eV is sufficient to destroy this matching 
even for ideal epitaxy. We checked this by calculating 
the conductance for a small bias voltage using the sur- 
face transmission function (STF) method introduced in 
Ref. 0- As expected, at 0.02 eV bias voltage the con- 
ductance becomes fully dominated by majority-spin elec- 
trons. Disorder would also tend to break the matching of 
the interface resonances even at zero bias. Therefore, we 
argue that in real Fe/MgO/Fe MTJs the minority-spin 
channel in the parallel configuration is closed. 

Unlike the parallel configuration, the interface reso- 
nances do contribute to the conductance in the antipar- 
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FIG. 2: Transmission probability as a function of kii . (a)- 
(c): Fe/MgO/Fe junction with 4 MgO MLs; (d)-(f): Same 
junction with Ag interlayers. (a), (d): Majority spins; (b), (e): 
Minority spins; (c), (e): Each spin channel in the antiparallel 
configuration. 




MgO thickness (ML) 

FIG. 3: Conductance (left axis) and TMR (right axis) 
vs barrier thickness for Fe/MgO/Fe junctions (open sym- 
bols). Triangles: majority spin, parallel configuration; dia- 
monds: minority spin, parallel configuration; squares: each 
spin, antiparallel configuration; circles: TMR ratio, calcu- 
lated disregarding minority spin in the parallel configura- 
tion (see text). Solid symbols: conductance and TMR for 
a Fe/Ag/MgO/Ag/Fe junction. 



allel configuration, where they tunnel into majority-spin 
states of the other electrode. The latter have no fine 
structure in the Brillouin zone, and hence the conduc- 
tance is weakly sensitive to a potential mismatch at 
the two interfaces which might occur in real junctions. 
The enhanced contribution of these interface resonances, 
which is clearly seen in Fig. ^p, leads to the decrease of 
TMR at low barrier thickness. We emphasize the fact 
that although the exact location of the interface reso- 
nances is not determined accurately due to intrinsic lim- 
itations of the den sity functional theory, their presence 
at the Fermi level [l3| inevitabl y re sults in the reduced 
TMR at small barrier thickness [lfj(. 

These features are evident in Fig. [3] which shows the 
conductance and TMR as a function of barrier thick- 
ness. In the parallel configuration the majority-spin con- 
ductance is controlled by the Ai band which dominates 
at large barrier thickness making TMR very large 0. 
Below 6 MLs of MgO, however, minority-spin electrons 
overcome the contribution from majority-spin electrons 
due to the contribution from the interface resonances. 
In the antiparallel configuration the spin conductance 
decreases faster than the majority-spin conductance in 
the parallel configuration, because it is dominated by the 
same interface resonances located away from the T point 
(see Fig. QJ. As was justified above, for real MTJs the 
minority-spin conductance in the parallel configuration 
can be disregarded in the calculation of TMR. This leads 
to the increase of TMR with increasing the barrier thick- 
ness. A similar behavior is observed experimentally until 
the barrier thickness exceeds approximately 1.5 nm [l(i| 
which corresponds to 7-8 MLs of MgO. At larger thick- 
ness the rate of decay for the parallel and antiparallel con- 
ductance becomes essentially identical. This crossover 
may be due to the loss of kn conservation induced by 



subbarrier scattering on defects, which makes tunneling 
electrons to diffuse over the surface Brillouin zone. The 
epitaxial junction model is inapplicable in this regime. 

In order to enhance TMR for thin MgO barriers we 
propose to use thin epitaxial Ag interlayers deposited 
at the Fe/MgO interfaces. Since the lattice parameter 
of Ag is close to both Fe and MgO lattice parameters, 
Ag can be deposited epitaxially on Fe(OOl) |25j, Fe can 
be grown on Ag j2?j], and Ag on MgO [23. Therefore, 
epitaxial Fe/Ag/MgO/Ag/Fe(001) tunnel junctions are 
feasible. It is known that an epitaxial Ag overlayer on 
Fe(OOl) surface notably modifies the electronic structure 
of the surface states 25], and it is natural to expect simi- 
lar changes for the Fe/MgO interfaces where Fe and MgO 
interact only weakly. If the minority-spin interface DOS 
is reduced by Ag, the antiparallel conductance will be 
suppressed. On the other hand, the majority-spin con- 
ductance should not strongly be affected due to almost 
perfect transmission through the Fe/Ag(001) interface 
[26j. This is the rationale for using Ag interlayers. 

We place 1 ML of Ag atoms on each Fe(OOl) electrode 
in the 4-fold hollow sites. The 4 ML MgO barrier is in- 
serted between Ag-terminated electrodes so that O atoms 
at the interfacial ML of MgO lie above the Ag atoms. 
This interface structure is considered the most stable for 
Fe/Ag(001) and Ag/MgO(001) interfaces 0, H|. To 
find the equilibrium interlayer distances, we relax the 
atomic structure of the MTJ using the pseudopotential 
plane-wave method [3(J implemented within the Vienna 
Ab Initio Simulation Package (VASP) [3l| . The general- 
ized gradient approximation 32] is used for the exchange- 
correlation energy. We find a 5.2% reduction in the Fe 
interlayer distance at the interface, the distance between 
the interface Fe and Ag layers being 1.88 A, and the dis- 
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tance between Ag and MgO layers being 2.76 A. 

Figs.|2Jl-f show the k||- and spin- resolved conductance 
of Fe/MgO/Fe junctions with Ag interlayers. Not un- 
expectedly, the majority-spin conductance is weakly af- 
fected by the Ag interlayers, whereas the minority-spin 
conductance and the spin conductance in the antiparallcl 
configuration change dramatically. The most pronounced 
difference for the latter two is the disappearance of the 
interface resonances that dominated the conductance of 
the Fe/MgO/Fe junction with no Ag interlayers (compare 
Figs. 12b, c and Figs.01,f). This strong change occurs due 
to the Fe-Ag hybridization which makes the interface res- 
onant band more dispersive and hence removes the Fermi 
level crossing responsible for the highly conductive reso- 
nant states. A careful examination of the band structure 
shows that the interface resonant band still crosses the 
Fermi level very close to the F point (an obscure circular 
feature in Figs.|5^,f), but due to its dispersive nature the 
interface DOS is small. As a result, this band crossing 
contributes 30% of the total minority-spin conductance 
in the parallel configuration, and only about 7% of the 
conductance in the antiparallel configuration. The sig- 
nificant reduction of the conductance in the antiparal- 
lel configuration leads to dramatic enhancement of TMR 
which changes from about 130% to 930% (see Fig. I3J. 
Thus, Ag interlayers practically eliminate the contribu- 
tion from the interface resonances and therefore enhance 
TMR for thin barriers. 

In conclusion, we have found that interface resonant 
states in Fe/MgO/Fe(001) tunnel junctions contribute 
to the conductance in the antiparallel configuration and 
are responsible for the decrease of TMR at small bar- 
rier thickness, which explains the experimental results of 
Yuasa et al. [lOj. Depositing thin Ag interlayers at the 
Fe/MgO interfaces is an efficient and practical way to 
suppress the tunneling conductance through these reso- 
nant states and thereby to enhance the TMR for thin 
barriers. 
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